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Abstract 

We report an experimental study on the dynamics of a thin film of 
polymer solution coating a vertical fiber. The liquid film has first a con- 
stant thickness and then undergoes the Rayleigh-Plateau instability which 
leads to the formation of sequences of drops, separated by a thin film, mov- 
ing down at a constant velocity. Different polymer solutions are used, i.e. 
xanthan solutions and polyacrylamide (PAAm) solutions. These solutions 
both exhibit shear-rate dependence of the viscosity, but for PAAm solu- 
tions, there are strong normal stresses in addition of the shear-thinning 
effect. We characterize experimentally and separately the effects of these 
two non-Newtonian properties on the flow on the fiber. Thus, in the flat 
film observed before the emergence of the drops, only shear-thinning ef- 
fect plays a role and tends to thin the film compared to the Newtonian 
case. The effect of the non-Newtonian rheology on the Rayleigh-Plateau 
instability is then investigated through the measurements of the growth 
rate and the wavelength of the instability. Results are in good agreement 
with linear stability analysis for a shear-thinning fluid. The effect of nor- 
mal stress can be taken into account by considering an effective surface 
tension which tends to decrease the growth rate of the instability. Finally, 
the dependence of the morphology of the drops with the normal stress 
is investigated and a simplified model including the normal stress within 
the lubrication approximation provides good quantitative results on the 
shape of the drops. 
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1 Introduction 

The fiber coating process is widespread in numerous industrial applications such 
as the manufacture of glass, polymeric and optical fibers, conducting cables or 
textile fibers. The application of a thin layer on these solid substrates should 
ensure mechanical or optical properties of the final deposit. Hence, it is of crucial 
interest to control the final thickness of the liquid film. It has been well known 
since [27] that a cylindrical free surface of a fluid is unstable under the action of 
the surface tension. Later, [1] described in his monograph the patterns observed 
in a spider web with a sticky fluid. He reproduced the experiment with castor 
oil and a quartz fiber, and detailed the spatial variations of the film thickness. 
For a fiber drawn out of a bath, [18] and later [30] provided a first view of the 
flowing regime as a function of the capillary number and the Goucher number 
defined as the ratio between the radius of the fiber and the capillary length. 
Subsequently, different configurations have been studied to identify the different 
mechanisms responsible for destabilization of the film and the dependence on 
various parameters such as the radius of the fiber, the viscosity, the inertial 
forces or the influence of surfactant on the growth rate of the instability [TH] H7J 
25, 26, 6 , Uni EJ [5H HI! ]■ More recently, we have reported a flow regime diagram 
which identifies, depending on the fiber radius and the flow rate on the fiber, the 
dominant physical mechanisms [lit fl~2 . 10, ]. It appears that for small fiber radii 
compared to the capillary length and low flow rates, the liquid film is dominated 
by the surface tension, and the instability mechanism is the Rayleigh-Plateau 
instability whose nature is absolute. For higher values of the parameters (fiber 
radius and/or flow rate) two other regimes have been discerned respectively 
dominated by gravity and drag (the drag/dravity regime) or by inertia and 
drag (the drag/inertia regime). The nature of the instability is then convective 
for these regimes. Extensive theoretical studies have investigated the dynamics 
of the film [151 [251 HOI ]■ Among the most recent studies, we can mention the 
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works of [21] and [5] on film thickness of the same order as the radius fiber 
and with negligible inertia contribution. Their numerical results are compared 
with experimental results ] performed in the three kinds of regime depicted 
above with a predictable deviation for the Drag/Inertia regime. 

In most industrial situations, the coating fluid is a polymer material or a 
complex fluid and exhibits non-Newtonian properties depending on miscella- 
neous parameters such as concentration, structure or flexibility of polymers. 
Yield stress, shear-thinning or clastic effects are some of the non-Newtonian 
behaviors which can affect the structure of the flow, the appearance of the in- 
stability or the morphology of the patterns. jS], in the case of a "dip-coating" 
configuration where the fiber is drawn out of a bath of liquid, observed that 
the film swells due to the presence of polymer in the solution. Considering the 
normal stress and the lubrication approximation, they found an analytical ex- 
pression of the film thickness as a function of the withdrawal velocity and the 
normal stress coefficient. In general for most of the cases, when the instability 
is studied for non-Newtonian fluids it reveals different classes of patterns and 
the presence of polymers can drastically change the dynamic of the system, as 
in Faraday or Saffman- Taylor instabilities 22> ]. In the case of the instability 
of a liquid jet [7J EH E ], the addition of polymers causes the formation 
of a "beads-on-a-string" structure where adjacent beads are joined by a thread 
which grows thinner and strongly delays the detachment of droplets. In this 
configuration, the flow is subject to a strong elongation: a velocity gradient 
exists in the direction of the flow due to gravity forces. This additional resis- 
tance to breakup compared to a simple fluid is due to large extensional stresses. 
Contrary to the case of a liquid jet, there is no elongational viscosity in the case 
of the flow down a fiber due to the no-slip condition on the fiber as detailed 
further. 

A complication inherent in the use of these complex fluids is that they ex- 
hibit different non-Newtonian properties with opposite effects. Notably, most 
polymer solutions are both shear-thinning and present elastic effects. In the con- 
text of fiber coating, and to obtain independently the role of the shear-thinning 
effect and normal stress on the Rayleigh-Plateau instability, we have performed 
experiments with two different polymer solutions: one with a rod-like poly- 
mer (xanthan), exhibiting a pure shear-thinning effect; the other with a flexible 
polymer (polyacrylamide abbreviated as PAAm) which exhibits non-negligible 
normal stress along with shear-rate dependence of viscosity similar to xanthan 
solutions. 

Our experiments are all performed in a regime where inertial and gravita- 
tional forces are negligible compared with capillary, elastic and viscous forces. 
The Rayleigh-Plateau instability is then absolute (T2J ] and the flow patterns 
consist of drops, where fluid is partly trapped in a recirculation zone, sliding 
down a very thin (smaller than 100 /im) and quasistatic liquid substrate. This 
droplike wave train emerges from a constant film thickness, i.e. the flat film 
region. We investigate the role of the non-Newtonian properties on such a flow, 
i.e, the flat film and the drop- like wave train resulting from the Rayleigh-Plateau 
instability. 
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This paper is organized as follows. In section 12. 1[ we proceed to a rheological 
characterization of the solutions. In section[2j we present the experimental setup 
and visualisation techniques. In section [31 we look at the flat film before the 
appearance of the instability and the shear-thinning effect on the film thickness. 
In section^ the experimental growth rate and the wavelength of the instability 
are measured experimentally for xanthan solutions (only shear-thinning effect) 
and for PAAm solutions exhibited strong normal forces and similar viscosity 
shear-rate dependence. These experimental data are then compared to the re- 
sults of a linear stability analysis taking into account a non constant viscosity 
with the shear-rate. In the last part, section[5j we will treat the effect of normal 
stress on the morphology of the drops and provide a simplified model to explain 
the dependence of the drop shape with the normal stress. 

2 Materials and characterizations 

Before considering the non-Newtonian properties, certain conditions are re- 
quired: first, to avoid inertial forces, and second to ensure a perfect wetting 
on the fiber. In a previous paper [TTJ ], we have presented a diagram of the 
expected flow regime which details the dominant physical mechanisms in the 
plane of the dimensionless numbers R/l c versus ho/R with R, the radius of the 
fiber, l c the capillary number and ho the flat film thickness. To ensure that the 
flow is dominated by capillary forces, with no inertia, some conditions on the 
fiber radius, the surface tension and the viscosity of the fluid must be fulfilled 
in agreement with the flow regime diagram mentioned above. In that capillary 
region, the flow on the fiber consists of drops sliding on a quasistatic thin film. 
To ensure a good wetting on the fiber and consequently an axisymmetric pat- 
tern on the fiber, the surface tension has to be lower than 40 mN/m. Finally, 
these two conditions are satisfied by using a fiber radius equal to 0.28 mm and 
solutions composed of a mixture of water, glycerol (to increase viscosity) and 
surfactant to reach a surface tension close to 30 mN/m. 

Some previous experiments carried out with Newtonian fluids on the same 
experiment [TTJ ], have shown that film thicknesses between 0.1 and 1 mm 
are possible with velocities ranging from 1 to 10 cm/s and viscosities between 
50 and 500 mPa.s. This implies that the shear-rate range encountered in our 
experiments is from 10 to 1000 s _1 . 

The liquids used are then semi-dilute solutions of polymers: xanthan and 
Polyacrylamide (PAAm) purchased from Sigma-Aldrich. Both present shear- 
thinning behavior but only the second one exhibits large elastic effects in the 
shear-rate range considered in the experiments. To discern the effect of elasticity, 
Boger fluids would have been theoretically more appropriate. Nevertheless, in 
practice, the large quantities of liquid required for our experiments and the 
possible degradation of the nylon fiber by the solvents used in Boger fluids, are 
two reasons prohibiting their use. 
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2.1 Samples preparation 



2.1.1 Polyacrylamide 

PAAm is a high molecular weight chain resulting in long flexible chains (M ~ 
1.5 x 10 7 g/mol). The samples are prepared very carefully to get homogeneous 
solutions. First, a solvent is prepared by mixing 50% of purified water and 50% 
of pure glycerol (all percentages presented in this article are weight percentages). 
Six polymer concentrations were studied in the range from 0.1% to 0.6%. The 
surfactant selected to reduce the surface tension is Triton X-100 (TX-100). This 
choice was motivated by its mixing properties at high polymer concentrations 
and for the resulting low surface tension [32, ]. The TX-100 concentration is 
4.5% (about 300 times the Critical Micelle Concentration, CMC, in water). 
For this high concentration we assume that the surfactant mobility timescale 
is higher than the timescale for the instability growth rate. Indeed, the time 
variation of the surface tension due to diffusion of surfactant is 10 -2 s for a 
TX-100 concentration 50 times the CMC in water [T3]. This time scale is of 
the same order of magnitude as the characteristic growth rate of the instability: 

^ ~ soxio-^o^io-a) 3 = °- 05 s ( usin § typical values for the viscosity 77, the 
fiber radius R, the surface tension 7, and the film thickness h). Consequently, 
in the range of concentrations considered in our experiments we assume that 
the interface is rapidly saturated with surfactant molecules before the instability 
occurs. The surface tension of the final solution is then 7 = 32.3±0.5 mN/m. We 
should note that no apparent rheological modifications are observed by varying 
the TX-100 concentration. 



2.1.2 xanthan 

xanthan is a rigid rod-like polymer (M ~ 5 x 10 6 g/mol). A preparation proto- 
col similar to the PAAm solutions was used for xanthan solutions. The solvent 
is slightly different: 60% of glycerol and 40% of water. Since xanthan is a 
polyelectrolyte polymer, the resulting rheological properties of this polymer are 
known to be modified by the addition of salt [311 ]• Thus, different concentra- 
tions of NaCl allow adjustment of the rheological properties. After the addition 
of TX-100, the surface tension of the final solution is 7 = 32.7 ± 0.8 mN/m 
independent of the salt concentration. 



2.2 Rheological characterization 

The rheology of polymer solutions was performed using an "Anton Paar" rheome- 
ter with a cone and plate geometry. We have chosen a large cone (radius: 49.988 
mm) with small angle (angle: 0.484°) in order to measure precisely normal 
stress in a large range of shear-rates. The sample temperature was fixed at 
20.00 ±0.05°C. 
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2.2.1 Polyacrylamide 

The evolution of the apparent shear viscosity, 77, is plotted versus the shear-rate, 
7, for two concentrations in figure[T] As usual, rj decreases with 7 and, at a given 
shear-rate, increases with the polymer concentration. Such curves are typical for 
shear-thinning fluids where a constant Newtonian low-shear viscosity is followed 
by a power-law dependence before reaching the viscosity of the solvent at high 
shear-rates. These measurements can be reasonably fitted by a four-parameter 
Carreau model [2H ]: 

i) = Voo H — (1) 

(l + H) 2 ) -1- 

where 770 and rjoo are respectively the viscosity for the zero-shear limit and 
infinite shear-rate, and r denotes a characteristic time scale that measures the 
scale at which the shear-thinning effect becomes important. The exponent n is 
the power of the following Ostwald power-law equation: 

v = PT- 1 (2) 

The zero-shear limit increases rapidly with the polymer concentration, that 
is typical of entangled polymer solutions (see the inset of figure [TJ . Recently, 
it has been shown |23[ ] that in a good solvent the entangled concentration c e 
for PAAm is about nine times the crossover concentration c* ~ 0.2 g/L. The 
temperature effect on the samples indicates that the viscosity decreases by 10% 
for a 5°C increase. 



The normal stress measurements are presented in figure 2(a) as a function 
of the shear-rate for different PAAm concentrations over a wide range of shear- 
rates. A significant increase of the normal stress with the shear-rate is observed 
in accordance with: 

Ni = Vi7 2 (3) 
where is the first normal stress coefficient characterizing the fluid [21 ]. 



Values are given in figure 2(a) for several Polyacrylamide concentrations. In 
the inset graph, we have presented the data using a double log-plot. It appears 
more clearly that, except for PAAm solutions at 0.6%, there is a discrepancy 
between the data and the curve fit for low shear-rates (lower than 200 s _1 ), 
indicating dependence of the first normal stress difference with the shear-rate. 
Thus the 0.6% solution would be the best candidate for studying the normal 
force effect in section [5] Normal stress magnitude can be compared to viscous 



stress by estimating the ratio as a function of the shear-rate 2(b) This 

ratio increases with the shear-rate, highlighting the importance of the normal 
stress, which starts to be dominant compared to the shear-thinning effect for 
shear-rates larger than 100 s _1 . One should then expect a large amount of 
normal stress in the drops, for which the shear-rate is always larger than 100 
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Figure 1: Variations of the shear viscosity 77, of PAAm solutions vs. shear-rate 
7, in a log-log scale. Each symbol refers to a different polymer concentration. 
Plain curves correspond to data fits with a Carreau model. The inset shows the 
zero shear-rate viscosity 7? , as a function of the concentration in PAAm. 
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Figure 2: (a) Log-log plot of the normal stress, iVi, as a function of the shear- 
rate, 7, for PAAm solutions in water-glycerol (50% : 50%) solvent; 4.5% TX-100 
surfactant was added to the solutions. Plain curves correspond to data fits using 
equation El (b) versus shear-rateg7 , for PAAm solutions 




Figure 3: Viscosity 77, as a function of shear-rate 7, for two solutions constituted 
by water-glycerol (50% : 50%), 4.5% TX-100 and 0.4% xanthan. Red cross 
correspond to free-salt solution and blue one are for a salt concentration up to 
0.8%. Numbers indicate the slopes of the Ostwald power-law model from fits 
over a 7 range: [10; 2000] s _1 . 



2.2.2 xanthan 

Rheological measurements are typical of shear-thinning fluids where a constant 
Newtonian low-shear viscosity is followed by a power-law dependence before 
reaching the viscosity of the solvent at high shear-rates (figure [3]). Nevertheless, 
as shown in figured the power-law behavior failed to fit the experimental results 
in the whole range of shear-rates. Since xanthan is a polyelectrolyte, the solution 
rheology and molecular configuration are greatly affected by the solution's ionic 
strength. Thus by adding 0.8% NaCl to the solution, the shear-thinning effect 
can be adjusted to n = 0.73 in a reasonable range of shear-rates from 5 to 
2000 s _1 . No significant normal stress has been detected for xanthan below 
4000 s _1 , a shear-rate which is not expected to be reached in the experiment. 
Subsequently, we will exclusively use salted xanthan solutions as pure shear- 
thinning solutions. 



9 




(a) (b) z ▼ (c) 

Figure 4: (a) Experimental setup. Scheme showing a fluid flowing down on a 
fiber from a upper tank. The flow rate is controlled by a valve and guided with a 
nozzle, (b) Picture snapped by a high-speed camera with a telecentric lens (1 x ) . 
The white bar length is 5 mm. (c) Spatiotemporal evolution of the film obtained 
by a vertical linear camera passing through top drops (black lines). Letters A, B 
and C denote respectively the fiat film, the ordered and the disordered pattern 
regions. 



2.3 Experimental setup 



As depicted in figure |4(a) the fluid flows from an upper reservoir (diameter: 



14 cm) down a nylon fiber (diameter: 0.56 mm). The relative pressure variation 
is about 0.001% during one minute for the highest measured flow rates. The 
flow rate is controlled by a valve composed of two axisymmetric cones. The 
mass flow rate Q is measured from the weight variation of a collecting tank 
recorded by a computer-controlled scale. A transparent nozzle guides the fluid 
on the fiber. Its verticality is crucial to obtaining an axisymmetric flow and it is 
ensured by a mechanical device which enables very accurate fiber displacements 
with a sensibility of 2.4 arc sec. Two perpendicular cameras with zoom lens 
help to control the axisymmetry of the film flowing down the fiber. 

As depicted in figure |4(b)| the flow presents three regions along the fiber. 
A meniscus is followed by a flat film with a constant thickness on a distance 
called the healing length, which increases slightly with the flow rate [HI ] . Then, 
the Rayleigh-Plateau instability leads to the formation of a regular pattern of 
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beads flowing on a very thin and flat film. In this paper, we will deal only 
with regime dominated by capillary forces, so we exclude low flow rates (the 
dripping regime) and high flow rates where inertial flow dominates. The flow 
regime is then absolute. The film thickness and the shape of the drops are 
captured by a high-speed digital camera with a telecentric lens. The interface 
position is detected in both space and time, so we are able to measure the 
film thickness h(z,t) with an accuracy of 0.02 mm. A linear camera provides 
spatiotemporal diagrams which deliver information on the dynamic of the flow. 
A vertical pixel line passing through the peaks of the drops is recorded and stored 
at constant time intervals. The resulting spatiotemporal diagram produces the 
(z, t) trajectories of the drops along the fiber. A typical spatiotemporal diagram 



is shown on figure 4(c) The uniform grey region, located at the upper part of 
the fiber is the place where the film is flat (region A) and gives rise to a zone of 
regular stripes with a constant wavelength and velocity for the drops (region B). 
Finally, downstream, some coalescences between drops lead to the formation of 
a disordered pattern (region C). 



3 Flat film 

In this section we focus on the region close to the inlet where the film thick- 



ness is constant (flat film, grey uniform region, see figure 4(a) I. In the case of 
Newtonian fluids, the thickness of the film, h, is given by the classical Nusselt 
solution [12l ]• I n the case of very thin films (h <C R), i.e. the planar case, 
there is a cubic relation between the flow rate on the fiber and h. We define 
the cylindrical coordinates system (r, 9, z), where r is the radial coordinate (the 
fiber center is the origin), 9 the azimuthal coordinate and z the axial coordinate 
oriented downward in the flow direction. In the case of a shear-thinning solu- 
tion exhibiting normal stress effects, for a steady axisymmetric flow, the stress 
balance in the axial direction z is written as: 

da zz ld(ra rz ) dp 

~5 1 a = "a P9 4 ) 

Oz r or oz 

for R < r < R + h{z), where a denotes the stress tensor, p is the pressure 
field in the film, p and g are respectively the fluid density and the gravitational 
acceleration. Since the first normal stress difference, N\, can be expressed as 
®zz — &rr = ipi (§7) i with v(r, z) the axial velocity, which varies along the film 
thickness, equation (Q} becomes 

diVi 1 d(ra rz ) dp da rr 

~a h « = « « P9 ( 5 ) 

oz r Or Oz Oz 

As no spatial variations are detected in the flat film, z-invariance of the 
velocity field implies that = 0. Thus, normal stress has no effect on the flat 
film. 
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To calculate the velocity profile in the flowing film, we model the shear- 
thinning effect by a power-law in accordance with expression [2j Thus ([5]) be- 
comes 

where the following dimensionless variables are introduced: f = h — 



l/n 



R ' H' 



The fluid velocity satisfies two boundary conditions: no-slip on the fiber 
(v(f — 0) = 0) and zero tangential stress at the liquid-air interface (drv(r = 
h) = 0). Considering this last boundary condition, equation (j6|) becomes: 



dv 
dr 



(j^(h(l + h/2)-f(l + f/2))y (7) 
The flow rate per unit length q = 2 ^ p r i s written in dimensionless form as 

1 f h 

q = q— = J v(r + l)df (8) 

We note that for the Newtonian case n = 1 and f3 ~ i], we recover the 
analytical Nusselt solution v = g| (2(R + h) 2 In (£) - (r 2 - R 2 )). 

Equation ([7]) is solved using a Runge-Kutta algorithm (starting at v(r = 0) = 
in order to satisfy the boundary condition on the fiber). First, the influence 
of the parameter n is studied for a constant flow rate q = 1, and the integral 
of the equation [8] is estimated by the trapezium rule. We choose two values 
for the h parameter (h 1U i n and h max ) satisfying q(h min ) < 1 < q{h ma _ x ) and we 
find the film thickness by a bisection method for which the condition q = 1 is 
satisfied at 0.1%. The results are shown in figure [S] The numerical solution 
for the Newtonian case, n = 1 is identical to the analytical Nusselt solution. 
For a constant flow rate, an increase in the shear-thinning effect modifies the 
velocity profile shape: the parabolic profile tends to be replaced by a plug- 
like profile. This results in higher velocity gradient close to the fiber, whereas 
close to the interface, the velocity gradient is almost zero. The film thickness 
is always smaller than for Newtonian fluids. The film thickness is plotted as a 
function of the flow rate on the film for a PAAm solution (0.4%, n=0.71), in 
figure ([6]). We choose to present experimental data only for PAAm solutions, 
since for xanthan solutions, the healing length is two or three times smaller than 
for PAAm solutions. The good agreement between the numerical solutions and 
our experimental data validates the choice of the Oswald power-law model for 
the viscosity and also our assumption on the negligible effect of surfactant on 
the zero-shear stress boundary condition at the liquid-air interface. 
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Figure 5: Velocity profiles of the flowing film at a constant flow rate (q = 1) for 
several n values from a Newtonian fluid (n = 1) to a high shear-thinning effect 
(n = 0.4). The inset is a close-up of the liquid-air interface region. 
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Figure 6: Numerical solution (solid line) and experimental results for a con- 
centration in PAAm equal to 0.4%. The dashed curve is the analytical solution 
for h <C R given by the equation H"6l for n = 0.73 . 



14 



4 Instability growth rate 

The impact of the shear-thinning and elastic effects on the growth rate of the 
Rayleigh-Plateau instability is investigated through experiments with xanthan 
(0.8% NaCl) and PAAm solutions (0.4%) in order to distinguish the role of 
each non-Newtonian property. After the flat film region, some variations on 
the film thickness are detected and a regular pattern of drops emerges due 



of the drop-like pattern is plotted in figure [7] as a function of the distance to 
the entrance nozzle. The wavelength increases (regime A in figure [7j until it 
reaches a well-defined value (regime B in figure [7]). Then, lower down, some 
coalescence events can disrupt the regular pattern (regime C particularly in 
figure [7^,). The wavelengths of the regular pattern for non-Newtonian fluids are 
somewhat higher than those expected with Newtonian fluids. Nevertheless, in 
both cases, the classical Rayleigh Plateau wavelength fails to fit the experimental 
data and is always smaller. The length of the regular pattern depends on the 
flow rate (at high flow rates, coalescence events occur earlier) but it is typically 
is of the order of seven centimeters for PAAm solutions (figure [7^,) and shorter, 
about four centimeters, for xanthan solutions (figure [TJd) . In the former case, 
we can note that the axisymmetric conformation is not the only case observed 
on the fiber: some non-axially symmetric conformations can be observed with 
asymmetric drops. Such conformations have been described by [5] as a roll-up 
transition and must be avoided in our case. 

In order to characterize the instability growth rate, we record a stack of 
images at a typical frame rate of 1000 images per second. Then, the position 
of the film interface is detected over space and time: h(z,t). Figure [5] shows 
the average film thickness over the time < h(z,t) > t and the extremal film 
positions for a PAAm solution. It shows successively the meniscus, the flat 
film and the onset of the instability which is marked by a strong variation of 
the film thickness. The velocity of the interface is calculated for each stack. 
Then, a point on the interface (chosen to become a point of maximum height) is 
followed at this velocity using the set of data h(z, t). The resulting values of the 
normalized profile (h — ho) /ho as a function of time for a typical experiment are 
plotted in the inset of figure [5] and fitted by an exponential law h ~ h '*" — Ae nt 
in the early linear stages. From this fit, the growth rate, fi, is extracted and 
averaged over several other experiments; the results for xanthan and PAAm 
solutions were reported in figure |9] 

A first simplified attempt to obtain an expression for the growth rate consists 
in a linear stability analysis. The fluid is assumed to exhibit pure shear-thinning 
effects with 77(7) = /37™" 1 , and we assume very thin films such that h <C R (pla- 
nar approximation). Thus, in cartesian coordinates, the following momentum 
equation holds, in the lubrication approximation: 




The wavelength 




(9) 



where II is the pressure gradient given by: 
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Figure 7: Wavelength at different flow rates for (a) xanthan and (b) PAAm 
solution. Letters A, B and C denote respectively the growth of the instability, 
the ordered and the disordered pattern regions. 
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Figure 8: Average film thickness in time < h > t along the z fiber axis (Q = 0.032 
g/s). Bars indicate the extreme values of the film thickness. The flat film 
thickness is denoted by ho. 
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Figure 9: The growth rate, ft, is plotted versus the flat film thickness, ho. The 
experimental data are represented by the dots and the equation [22] by curves. 
The inset shows the growth of the film using the method described in section 21 
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n = pg + 7 



d z h ^d 3 h\ (d z h d 3 h\ 



Given that there is no fluid slippage on the fiber and no stress on the liquid- 
air interface, the velocity is calculated from the momentum equation, so that 



v(r) 



1 



n 



l/n 



h W/n _( h _ r) l+l/r, 



1 + l/n V/3, 

The flow rate per unit length, defined as q p — f£ vdr, is given by: 

/n\ 1/11 fe 2+1 /" 

qp -{f3) 2 + l/n 
and also satisfies the mass conservation equation: 

dh , d % = 
dt dz 



(11) 



(12) 



(13) 



We assume infinitesimal perturbations around the uniform film thickness, 
ho, and for the corresponding flow rate, q p o, so that 



h(z,t) = h + hi(z,t), 
% = <7p0 + %i, 

This results in the following expressions 



(14) 
(15) 



<Z P o 

q P i 



, l/n ,2+1/r. 

P9\ K 



/3 J 2 + l/n 

7 / d z hi d 3 hi 



npg \ R 2 

and in the linearized equation: 

_ dhi _ <7 p o7 ( d 2 z h x t d 4 h x 
dt 



dz 3 



2 + - IT 

n J h 



2 + 



1 \ q P o dh! 



n J ho dz 



(16) 

(17) 
(18) 



(19) 



npg \ R 2 dz 4 

Developing the thickness perturbation as hi = Ae llykz ^^ leads to the dis- 
persion relation: 



/i V n / n P5 V# 



The maximum of Im(w(fc)) gives: 



9po7 
4npgR 4 '' 



h<^R 



(20) 



(21) 
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As shown in figure [51 the flow rate given by dadoes not reproduce the exper- 
iment where the planar approximation (h <C R) is not valid. In previous works 
on Newtonian fluids, it has been shown [5] ] that for h ~ R, the expression for 
the growth rate is similar except that R should be replaced by R + h. More- 
over, we made the choice to use our numerical calculation described in section[3] 
which provides g num (/i) without assumption on the film thickness. Finally, the 
addition of a large amount of surfactant modifies the growth rate by a factor 4, 
as described by [5] since surfactants change the surface elasticity. 

Taking into account these corrections, we obtain the expression 

1 9nu m (fe)7 (22) 

which is plotted in figure [9] for both chemical systems. 

Data for the xanthan solution are well described by [22] To validate our 
method and to compare with a Newtonian fluid of similar surface tension (7 = 
20.9 mN/m) and viscosity (77 = 0.965 Pa.s), we have performed an experiment 
with a silicon oil (with p — 96.5 kg.m~ 3 and ho = 0.55 mm) and measured a 
growth rate f2 equal to 10.9±0.7 s _1 . Equation [22] for n = 1 gives fi equal to 
11.1 s _1 . The small variations in growth rate between the xanthan solution and 
the silicon oil are reasonable as the value of the viscosity is of the same order. 
Concerning PAAm solutions, there is significant deviation from the theory due 
to the normal stress of this solution. A qualitative explanation of the role 
of normal stress can be provided by the "hoop stress" effect pjjj ]. A short 
description of this effect can be made by considering the liquid surface as an 
infinite cylindrical shell of thickness e ( figure [T0|) . For a cylinder of radius R + h, 
the balance between the internal pressure P and the stretching stress age, leads 
to 2(R + h)LP = 2Laeee. The internal pressure P is generated by the normal 
stress N\ in the bulk. Interpreting the stretching force per unit length, 2<rgge, 
in term of surface tension 7^ , we obtain 

7Vi = -ipij 2 (R + h) (23) 

An estimation for 7 = 100 s _1 gives a 7^ of about —10 mN/m significantly 
lowering the effective surface tension. The growth rate should therefore be 
estimated with the effective surface tension lower than the fluid surface tension 
and the resulting curve for the growth rate would be shifted and enable us to 
recover the experimental data. 

5 Drop morphologies: normal stress effect 

This section is devoted to a comparison between the patterns of flowing films 
of PAAm and xanthan, focusing particularly on the normal stress effect on the 
shape of the drops. Such a comparison requires polymeric solutions having 
similar shear-thinning properties. Further, 0.8% NaCl was added to xanthan 
(0.4%) solutions to decrease the high shear-thinning effect (figure [3]). Optimal 
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R + h 



Figure 10: Stretching stress, crgg, in a thin cylindrical shell of radius R + h and 
thickness e <C R + h. 



adjustment of the shear-thinning of PAAm solutions (0.6%) is achieved as shown 
in the inset of figure [TO Thus, the difference between the two solutions concerns 
only the presence or absence of normal stress. 

The typical pattern observed on the fiber consists of an axisymmetric film 
of constant thickness. Then the film breaks up spontaneously into a drop-like 
wave train as described in the previous section. 

For axisymmetric patterns, the superposition of drop profiles in a PAAm 
and xanthan films is shown in figure [TTJ Clear differences can be noticed in 
the profiles, notably the steepening of the drop front for the PAAm solution 
compared to the xanthan solution. In both profiles, there is a clear asymmetry 
between the front and back of the drops which is more accentuated for the 
xanthan drop. This remark suggests that the shape of the drops is affected by 
gravity. The apex heights of both drops are identical as well as the film substrate 
between drops (the trailing edge for PAAm is longer than for xanthan). These 
observations confirm the fact that the normal stress plays a significant role in 
the thin regions, close to the tail and the front of the drops, but exhibits no 
effect in the center of the drop (the thick region). To quantify experimentally 
the swelling effect observed with PAAm solution, we define the slope of the front 
H / L as shown in figure 1111 This parameter is plotted in PAAm and xanthan 
solutions for different flowing rates in figure [T2l 

To highlight the swelling process, we consider a scaling law analysis, starting 
from the stress balance [5] Since the film is not flat the z-invariance is no longer 
valid. Exhibiting the contribution of the normal stress difference a zz — a rr , and 
shear stress a rz , equation [5] becomes 

dQMtO 2 ) , l d(r V (m ^dp da rr 
dz r dr dz dz 

The axial velocity v(r, z), is determined using a series of the form [3J ]: 
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Figure 11: Drop shapes for two polymer solutions. The drop front of the vis- 
coelastic solution (dashed green line) is swollen compared with the pure shear- 
thinning liquid (solid red line). 
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Figure 12: Slopes of drop fronts for xanthan and PAAm as a function of the 
capillary number Ca. 
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v(r, z) — ao(z) + a\(z)r + a 2 (z)r 2 (25) 

Functions <Xq(z), a%(z) et 03(2) are calculated using the following three equa- 
tions, two for the boundary conditions and the last one for the mass conserva- 
tion: 

• Boundary condition at the interface with the fiber v(r = R, z) = 0; 

• Boundary condition at the liquid/air interface d^O— fl+ft( z )) _ q. 

• The mass conservation equation, 



dh d r R+h{z ' t] 

- + -/ „(r,*)dr = (26) 



R 



Considering equation[26]in the reference frame of a drop moving at a velocity 
and using the 

axial velocity is 



U and using the condition that the mean flow rate satisfies q ~ Uh — > 0, the 

h— >o 



3 (R-r)(r-2h( z)-R) 
2 2h 2 (z) 

Assuming a constant viscosity, 77, the stress balance equation[24]to zero-order 



m r is 



h(z) - R 2 9h' dp da rr 

ir>u ^mz) ~ u ^^(z-) = d- z ~^z-- p9 (28) 

The normal stress balance at the free surface of the him assumes that — P + 
a rr = "/k with k the curvature of the interface. If L is the characteristic length 
in the axial direction and if H is the characteristic apex height of the drop, 
then the normal stress balance at the free surface of the him, accounting for the 
curvature, is given by 

d P da rr _ J9g + ^h) (29) 



dz dz \ R 2 dz 3 

Thus the right-hand side of equation [28] becomes 

dp da rr ( H H , 2 



dz dz M~^[z? + Lff +l ') (30) 
where l c = is the capillary length. 

Experimental observations suggest that -p- -C jjp and k 2 -C jj^?- So, the 
scaling law analysis leads to the following equation: 

H 3 R 2 „ ( & 



T3~^al-6 T | (31) 
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where Ca = 2— is the capillary number and S£ = is the normal stress 
characteristic length. 

This scaling law gives the slope of the front of a drop, H/L, as a function 
of the viscoelastic properties of the polymeric solution. In particular, in the 
case of polymeric solutions exhibiting normal stress, ip\ ^ 0, the expression 1311 
clearly shows that H/L decreases. A comparison between the experimental data 
and the results of the scaling analysis is presented in figure [12] for different flow 
rates. There is a good agreement between the experiment and the model which 
succeeds in highlighting the swelling effect on the drop shape induced by the 
normal stress effect. 

6 Conclusion 

The effects of non-Newtonian properties of fluids have been investigated in the 
case of a film flowing down a vertical fiber. The flow on the fiber can be divided 
into three regions: (A) at the inlet, the film exhibits a uniform thickness, i.e 
the flat film region; (B) the uniform film is progressively replaced by a well- 
defined pattern of drops separated by a thin film, i.e the Rayleigh-Plateau region; 
(C) the coalescence of drops disrupts the flow and give rise to a disordered 
pattern. In order to disentangle the role of the shear-thinning effect and of the 
normal stress, we have considered two kinds of polymer solutions. The first 
consists of rigid rod- like polymers (xanthan), exhibiting a strong shear-thinning 
behaviour but negligible elastic effects. For the second solution, we used flexible 
polymers (PAAm) exhibiting strong elastic effects and shear-thinning effects 
similar to those of xanthan under certain physico-chemical conditions. Some 
adjustments have been made by modifying both the polymer concentration and 
the physico-chemical properties of the solutions to enhance or reduce one of the 
non-Newtonian properties : shear-thinning or elastic effect. 

Consequences of both effects have been investigated in the first two regions 
of the flowing film. In the flat film region, due to the invariance of the film 
thickness in the axial direction, only the shear-thinning effect is effective. At a 
constant flow rate, our experiments demonstrate that, as a consequence of the 
shear-thinning effect, the thickness of the film is always smaller than in the case 
of Newtonian fluid. Our results clearly show the influence of the shear-thinning 
effect on the velocity profile: a parabolic profile in the Newtonian case tends to 
become a plug-like profile. Thus, an increase of the shear-thinning effect yields a 
thinner, unperturbed film. Further downstream on the fiber, the film undergoes 
the Rayleigh-Plateau instability. The growth rate of the instability has been 
investigated experimentally and theoretically using a linear stability analysis. 
Good agreement is found between the experimental data for xanthan and the 
model. For PAAm solutions, and to take into account the normal forces, we 
consider an effective surface tension (lower than the fluid surface tension) which 
tends to decrease the growth rate and to recover the experimental data. The 
morphology of the patterns resulting from the instability depends on the non- 
Newtonian properties. In particular, the drops formed with PAAm solutions 
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exhibit a swelling effect compared to drops observed with xanthan solution, for 
a similar shear-thinning effect. We observe that the drop of fluid with normal 
forces is less rounded compared with the case of a pure shear-thinning drop. 
This swelling effect has been quantified by a scaling law analysis where the 
slope of the drop front is expressed as a function of the normal stress. 

As a conclusion, by considering two kinds of polymeric solutions with the 
same shear-thinning effect, which differ from each other in the presence of nor- 
mal forces, we have succeeded in understanding the relationship between the 
rheological properties and the destabilization of the flowing film on a fiber as 
well as the morphology of the observed patterns. This should be helpful in un- 
derstanding what happens with more complex fluids, in particular fluids which 
exhibit more elastic effects where the elasticity could prevent the growth of the 
instability. 
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